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Abstract. In our previous work, we have established the existence of transonic characteristic 
discontinuities separating supersonic flows from a static gas in two-dimensional steady compressible 
Euler flows under a perturbation with small total variation of the incoming supersonic flow over 
a solid right-wedge. It is a free boundary problem in Eulerian coordinates and, across the free 
boundary (characteristic discontinuity), the Euler equations are of elliptic-hyperbolic composite- 
mixed type. In this paper, we further prove that such a transonic characteristic discontinuity 
solution is unique and //-stable with respect to the small perturbation of the incoming supersonic 
flow in Lagrangian coordinates. 



CO 



1. Introduction 

We are concerned with the well-posedness of transonic characteristic discontinuities that separate 
supersonic flows from a static gas in two-dimensional steady compressible Euler flows. The govern- 
ing equations are the following Euler system that consists of conservation laws of mass, momentum, 
and energy: 

d x (pu) + dy(pv) = 0, (1.1) 
d x {pu 2 + p) + d y (puv) = 0, (1.2) 
d x {puv) + d y {pv 2 + p) = 0, (1.3) 
d x (puE + pu) + d y (pvE +pv) = 0, (1.4) 

where E = \{u 2 + v 2 ) + e. The unknowns p, p, e, and (u,v) represent the density, pressure, internal 
energy, and velocity of the fluid, respectively. Specifically, for a polytropic gas, the constitutive 
relations are 

p = p 7 exp(— ), e = —!—-, (1.5) 
c v 7 - 1 p 

where S is the entropy, c u is a positive constant, and 7 > 1 is the adiabatic exponent. The sonic 
speed c is determined by 

/ rln / r\tm 

(1.6) 
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In our previous work [5], we have established the existence of a weak entropy solution to the 
following initial-boundary value problem: 

CUE]) - Cd, in x > 0, y > g(x), 

< U = U , on x = 0,y > 0, (1.7) 

P = P, % = g'(x), on x > 0,y = g{x), 

provided that the incoming supersonic flow Uq is close in BV to a reference state U + , where 
U = (u,v,p,p), and U + = (u,0,p,p + ) is a uniform supersonic flow; that is, it is a constant vector 
and u > c+ = \J~^~- The unknowns in problem (jl.7p are the supersonic flow U = U(x,y) and 
the free boundary T>jj = {y = g(x), x > 0}, which is a Lipschitz curve passing through the origin. 
The free boundary is actually a transonic characteristic discontinuity (vortex sheet and/or entropy 
wave) separating the supersonic flow U from the unperturbed static gas U- = (0,0,p,p_) that is 
subsonic (cf. Figure [1]). Thus, problem (|1.7p is a free boundary problem (in the Euler coordinates) 
and, across the free boundary, the Euler equations are of elliptic-hyperbolic composite-mixed type. 



Solid Right- Wedge 



sonic Flow (I): (u,v,p,p) 




Vortex Sheet/Entropy Wave 



Static Gas (II): (0, 0,p,p. 



FlGURE 1.1. For supersonic flow passing the corner O with a static gas U- on 
the right of the solid right wedge (i.e., with velocity zero), then a combined vortex 
sheet /entropy wave is generated to separate the static gas below from the supersonic 
flow above. 

To further study the L 1 stability, we run into a difficulty to define the L 1 -distance between two 
solutions, since two solutions U and V describing the supersonic flows are generically defined in 
different domains depending on the location of their respective free boundaries T>u and T>y. For 
some given "time" x* > and at some point y* £ R, it may happen that, at (x*,y*), it is a static 
state for one solution, but is a supersonic flow for the other. To this end, we employ the special 
structure of the Euler equations of the two-dimensional steady flows, and introduce the Lagrangian 
coordinates (£,77) associated to a solution, as done in the study of transonic shocks [HI [Til [16], t° 
transform the free boundary to the positive £-axis, so that the solutions U and V are then defined in 
the same domain {£ > 0, rj > 0} in Lagrangian coordinates. Then the initial data functions Uq(0, y) 
and Vb(0, y) are transformed to the corresponding initial data functions C/o(0, 77) and Vq(0, rj) in their 
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respective Lagrange coordinates. Similar to the analysis in [5] for the construction of approximate 
solutions, we construct the approximate front tracking solutions U (£, 77) and V(£, 77) in the domain 
{£ > 0, rj > 0} with the initial data Uq{0, rj) and Vo(0, r/) on £ = 0, and the boundary data p = p on 
rj = respectively, and study the L 1 -stability property of these solutions for the initial-boundary 
value problem (|2.17p in Lagrangian coordinates, as proposed accurately below. With this merit 
of Lagrangian coordinates, we can show the characteristic discontinuity solutions are both unique 
for the given incoming supersonic flow in the Euler coordinates and L 1 -stable with respect to the 
small perturbations of the incoming supersonic flows in the Lagrangian coordinates. 

In Section 2, we present the Lagrangian coordinates associated to a solution to problem (jl.7p 
and formulate this problem in the Lagrangian coordinates as problem (|2.17|) . Then we briefly re- 
view how the existence of a weak entropy solution is shown by using the front tracking method. 
Section 3 is devoted to establishing the L 1 -stability of problem (|2.17p . Then, in Section 4, we show 
the existence of a semigroup associated with problem (|2.17p . Finally, in Section 5, we show the 
uniqueness of solutions to problem (|2.17p in the larger class of viscosity solutions. By Wagner's 
Theorem [151 Theorem 2], there is a one-to-one correspondence between a bounded measurable 
weak solution of a hyperbolic system of conservation laws and a bounded measurable weak solution 
of the corresponding equations in Lagrangian coordinates. Thus, the uniqueness results established 
for problem (|2.17p imply the uniqueness of transonic characteristic discontinuities in Eulerian co- 
ordinates. This is clarified in Section 6. The well-posedness results are summarized in Theorem 
16.21 These results on the existence, uniqueness, and L 1 -stability in Lagrangian coordinates, to- 
gether with the existence theorem established in [5], may be considered as a complete mathematical 
well-posedness theory of such transonic characteristic discontinuity solutions. 

2. Formulation of the Problem in Lagrangian Coordinates 

In this section we reformulate the problem in Eulerian coordinates to the problem in Lagrangian 
coordinates and show the existence of entropy solutions for the problem in Lagrangian coordinates. 

2.1. Lagrangian coordinates. For a piecewise C 1 smooth flow, we may introduce the Lagrangian 
transformation as used in [91 [16] to formulate problem (|1.7p in Lagrangian coordinates, which 
enables us to straighten the streamlines and hence treat a strict hyperbolic system derived from 

Let U £ L°°(M + ) and g E Lip(IR + ) be a piecewise C 1 weak entropy solution to problem (jl.7p . 
Define 

r(x,y) 

t] = rj(x,y;x ,y ) = pu(s,t)dt - pv(s,t)ds, (2.1) 

J (xo,yo) 

where (xo,2/o) is a fixed point on the transonic characteristic discontinuity T>u, and the integration 
is on any smooth curve T connecting (xo>2/o) with (x,y) and lies in the upper side of Vjj. Since 
the domain {x > 0,y > g(x)} is simply connected, by the conservation of mass, as well as the 
Rankine-Hugoniot (R-H) conditions for crossing the shock- front, 77 is a well-defined function of 
(x,y) with y > g(x), independent of the choice of T. Clearly, we have 
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We also note that, by the last condition in (|1.7|) . rj is independent of (xo>2/o) on Vjj. In fact, for 
(^0)2/o) and (x' ,y' ) on the characteristic discontinuity, we have 

f(x'o>y'o) 

rj(x,y;x ,y ) - i](x,y;x' ,y' ) = / pu(s, t) dt - pv(s, t) ds 

(pu{s, g(s)) ■ g'(s) - pv(s, g{s))) ds 



Xo 

= 0. (2.3) 

Since the characteristic discontinuity g(s) is in Lip(M + ,R), it is differentiable almost everywhere, 
so that the integrand is discontinuous only at these points in a set of Lebesgue measure zero, which 
is harmless for the above calculation. Hence, in the following, we may write r\ = r](x, y) with r/ = 
whenever y = g(x). 

If u, v, and p belong to L°°, and (jl.ip is satisfied in the sense of distributions, we can also find a 
unique Lipschitz continuous function i] so that r](xo,yo) = and (|2.2p holds. As a matter of fact, 
using standard mollification, we may approximate pu and pv in weak convergence of distributions by 
C°° functions (pu) € and (pv) e , for which the equality d x (pu) e + d y (pv) e = still holds. Furthermore, 
by the Young inequality, we have || (pu) e \\ Lac < C\\pu\\ LX > and IK/w^H^oo < C \\pv\\ LOO ■ Then, as 
above, we can solve if £ C 1 so that r] £ (xo,yo) = and hence, in any bounded domain, {rf} is a 
family of C 1 functions that are uniformly bounded and equicontinuous. Thus, there is a subsequence 
{r) k } that converges uniformly to some rj; by taking a diagonal subsequence, we can find r\ that 
is defined in the whole domain and r] k converges uniformly to ry in any compact subregion, hence 
r] k — > 7] in Qi 1 . It is obvious that r]{xo,yo) = 0, and (|2.2[) holds by uniqueness of limits in the sense 
of distributions. This then ensures that 7] is Lipschitz continuous. The uniqueness follows from 
the well-known fact that a distribution with zero derivatives must be a constant. Differentiating r/ 
along a streamline (which is Lipschitz continuous) and using (|2.2p yield that it is constant along 
the stream line. 

Now we introduce the following Lagrangian transformation (x,y) h-> (£, 77): 

(£,??) = (x, V (x,y)). (2.4) 

Then we have 




(2.5) 



d(x,y) 
thus 

d x = - pvd v , d y = pudrj. (2.6) 

This transform is Lipschitz continuous and one-to-one provided pu > 0. A simple computation 
shows that equations (|l.ip - (|1.3p may be written in divergence form: 

'*(-)- <*>(-) =0. 

\pu/ u J 

<%(« + £)-4i(£)=o, ^ 

d^v + d v p = 0, 
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or, as a symmetric system for U = (u,v,p) T , 

AdsrU + Bd v U = 0, 

with 
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(2.8) 



(2.9) 



For pu ^ 0, the conservation law of energy becomes d^^-^- — h ^rj) = 0, that is, 

r2 



u 2 + v 2 



+ 



b(r,). 



(2.10) 



2 ' 7-1 

This is the well-known Bernoulli law. As b(rj) is given by the initial data, in the following we focus 



on system (|2.T[) with p determined by U = (u,v,p) through (|2.10p . 
The eigenvalues A of (JZS]) (i.e. \XA - B\ = 0) are 



Ai 

A 2 
A3 



pc U (V 



v? — c 2 \u 



yM 2 



pc U (V 



u 2 — c 2 \u 



- + Vm 2 -i), 



(2.11) 
(2.12) 

(2.13) 



where M 



\f u 2 +v' 2 



is the Mach number of the flow. Then, for u > c, system (J2T8J) is strictly 



hyperbolic. The associated right-eigenvectors are 



n = + v, —it, -Aiii) T , 

P 

r 2 = (u,f,0) T , 

^3 = K 3 (— +V,-U,-\ 3 U) T , 
P 



(2.14) 
(2.15) 
(2.16) 



where Kj can be chosen that rj ■ VAj = 1, since the jth-characteristics fields, j = 1, 3, are genuinely 
nonlinear. Note that the second characteristic field is always linearly degenerate: r 2 • VA 2 = 0. 

2.2. Formulation of the problem in Lagrangian coordinates. As we noted that the transonic 
characteristic discontinuity becomes the positive £-axis in Lagrangian coordinates, we formulate the 
problem in Eulerian coordinates into the following initial-boundary value problem for equations 



2.17) 



(I2TD in f > 0,rj > 0, 

U(p,rj) = U (v) on e = 0,r ? >0, 
p = p on £ > 0, 7] = 0. 

Once we solved ?7 from this problem, we then obtain the free boundary in Eulerian coordinates 



g{x) 



u 



-(£,0)d£. 



(2.18) 



These are the corresponding forms in Lagrangian coordinates for problem (| 1 . 7j) . 
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2.3. Existence of entropy solutions. The entropy solutions of problem (|2.17p can be defined 
in the standard way via integration by parts. We note that the existence of entropy solutions to 
problem (|2.17p can be constructed easily by using the front tracking method as carried out in [5] 
provided that \\Uq — ^+Hbv(R+) * s sufficiently small. In particular, the following lateral Riemann 
problem with the boundary data p = p on the characteristic boundary {ij = 0} is uniquely solvable. 

Lemma 2.1. Consider the following lateral Riemann problem: 

JZH) in £ > 0, i] > 0, 
< U = U+ on ^ = 0,t)>0, 
p = p on £ > 0,7] = 0. 

There exists e > so that, ifU+ lies in the ball O e (U+) with center and radius e, then there is 
a unique admissible solution that contains only a 3-wave. 

Proof. 1. Note that system (12. 7ft is strictly hyperbolic for u > c. For each point U with u > c, in 
its small neighborhood, we can obtain C 2 -wave curves $>j(a;U),j = 1,2,3, so that 3>j(0;J7) = U 
and ^j^| a =o = rj(U): $j(at; U) is connected to U from the upper side by a simple wave of j'-family 
with strength \a\; for a > and j = 1,3, this wave is a rarefaction wave, while, for a < and 
j = 1,3, this wave is a shock. For j = 2, the wave is always a characteristic discontinuity. 

For our purpose, we note that there is also a C 2 -curve ^3(/3; U) which consists of those states 
that can be connected to U from the lower side by a 3-wave of strength (3. We have ^(O; U) = U 
and lI#l/3=0 = ~ r 3(U). 

2. We set U = (u,v,p) T in Lagrangian coordinates and use U3 to represent the third argument 
of the vector U. Then, to solve the lateral Riemann problem, it suffices to show that there exists 
a unique f3 so that (^3 (/?; U^))^ = p. Therefore, we consider the following function: 

L((3- u+) = * 3 (/3; U + )) 3 - vp 3 (0; U+)) 3 . 

It is clear that L(0; U + ) = 0, and 

dL(0;U+) - . - 

y Qp + ' = -rs(U + )3 = k 3 X 3 u + 0. 

By the implicit function theorem, there exists e > such that, for [/+ € O e (C/+), there is a 
function (3 = (3(U+) so that L(/3(t7+); U+) = 0. Then, similar to Lemma 2.2 in [5], by using the 
Taylor expansion up to second order (recall that is C 2 ), we obtain the following estimate: 

P = K\p + -p\ + 0{l)\U + -U + \\ (2.19) 

with a suitable constant K depending only on [/+. □ 

For the reflection of weak 1-waves off the characteristic boundary {rj = 0}, we have 

Lemma 2.2. Suppose that U l , U m , and U r are three states in O e (U + ) for sufficiently small e, with 
U m = $i(ai ;U l ) = $3(03; U r ). Then 

a 3 = -K 2 at + M 2 \a 1 \ 2 , (2.20) 
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with the constant K2 > and the quantity M2 bounded in O e (U_ + ). Furthermore, for U = 
(ui,vi,pi), \K%\ > 1, l-Ka I < 1, and \K?\ = 1 when vi < 0, v\ > 0, and vi = 0, respectively. 

The proof here follows the similar argument for Lemma 2.4 in [5]. 

For our current problem, we have to modify slightly the Glimm functional as done in [5]. In 
particular, this will allow us to show later that the Lyapunov functional $ (see Section 3) decreases 
when a weak wave of 1-family is reflected off the boundary. We now introduce the following version 
of the Glimm functional: 

0(O = v(£) + kG(O, 

where k > is a large constant to be chosen. The terms V and Q are explained below: 

• The weighted strength term V(£). We define the total (weighted) strength of weak waves in 
U s (£,-) as 

V(0=Z>«l- (2-21) 

a 

Here, for a weak wave a of i a -family, we define its weighted strength as 

J k + a if i a = 1, 
b a = < 

la if i a = 2, 3, 

where k+ > K2 is a positive constant which has been fixed by considering the reflection coefficient 
K2 in the strength 03 of the reflected 3-waves in the interaction between weak 1-wave and the 
characteristic boundary {r/ = 0} as given in Lemma 12.21 

• The interaction potential term Q(£). The interaction potential term we used here is a modifi- 
cation to the one introduced by Glimm, that is, 

2(0 = 1^1 + E IM = Q -4 + Qb, (2.22) 

where the set is of all the couples (b a ,bp) of approaching wave-fronts. The term Qt> in our 

wave interaction potential is an additional term, in comparison with the Cauchy problem; and if a 
weak wave a of 1-family is approaching the boundary, we write a G £^b- 

Now, for given initial data Uq, we may construct an approximate solution II s by a front tracking 
algorithm introduced by [12], where S is a small parameter measuring the accuracy of the solution, 
which controls the following four types of errors generated by the algorithm: 

• Errors in the approximation of initial data; 

• Errors in the speeds of shock, characteristic discontinuities (vortex sheet and entropy wave), 
and rarefaction fronts; 

• Errors by approximating the rarefaction waves by piecewise constant rarefaction fronts; 

• Errors from removing all the fronts with generation higher than N £ r L + (N depends on 5). 

The construction of a Glimm functional as above provides the necessary uniform estimates that 
guarantee the existence of a subsequence of II s which converges to a bounded entropy solution of 
flUTD in C([0,T];L 1 (M + )) for any T > 0. 
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The following sections are devoted to the L 1 -stability and uniqueness properties for problem 

(mm. 

3. A LYAPUNOV FUNCTIONAL AND L 1 STABILITY 

This section is devoted to establishing the L 1 -stability of two entropy solutions U and V of 
problem (|2.17p with respective initial data Uq and Vq obtained by the front tracking method: 

\\U{Z, ■) - V(£, -)|| L i (K+) <C\\U - Vb|| L i (M+) , (3.1) 

with a constant C depending only on the equations and the reference state Z7+. To this end, 
following an idea in [5J, we introduce a Lyapunov functional $>(U Sl , V s " 2 ) by incorporating additional 
new waves generated from the weak wave interactions with the characteristic boundary {tj = 0}, 
which is equivalent to the L-^R" 1- )— distance: 

Cr 1 \\U(t, •) - V(£, -)\\ L i < m V) < Ci||C/(^, •) - V(£, -)IUx (3-2) 

and decreases as £ increases: 

^H6r),^ 2 (6r))-^ 1 (6r),^ 2 (6,0)<C' 2 max(5 1 ,5 2 )(6-ei), V 6 > 6 > 0, (3.3) 

for some constants Cj, i = 1, 2, where and y^ 2 are two approximate solutions corresponding to 
the initial data Uq and Vo obtained by the wave- front tracking, with accuracy 5\ and 62, respectively. 

Estimate (|3.3p has many important implications. Firstly, taking Uq = Vq and £1 = in (|3.3p 
and using (|3.2p . we obtain lirn,5 1)1 5 2 _ ) .o H^ 1 ^ ') ~~ U S2 (£, ')\\ L im+\ = 0- Thus, the whole sequence 
of approximate solutions {U s } converges to the same limit. Secondly, taking directly 61,62 ~~ ^ 
yields (|3.ip . This again implies that the entropy solution of problem (|2,17p obtained by the front 
tracking method is unique. 

In the following, we first define specifically the Lyapunov functional and then verify (j3.2|) and 
Q. 

3.1. Definition of the Lyapunov functional and its equivalence to the L 1 distance. Simi- 
lar to [Hll3tH3]. when £ is fixed, for each 77 > 0, we solve a "Riemann problem" with the below-state 
f7(£, 77), the upper-state V(^,rj), by moving along the Hugoniot curves Si, C2, and S3 of system 
()2.7p in the phase space. We use hi(£,r]) to denote the strength of the i-Hugoniot wave in this 
"Riemann solver" , which is totally determined by U(£, rj) and V(£, rj). As in the standard Riemann 
solver, Yli=i is equivalent to \U(£,rj) — V(£,r])\, with a constant depending only on U + . 

Now we define the weighted L 1 -strengths: 

<k(t,v) = <$hi(M, (3-4) 

where the constants cf are to be determined later, depending on the estimates on the wave inter- 
actions and reflections off the characteristic boundary {77 = 0} (as done in Section 2). 
Then we define the Lyapunov functional: 

3 00 

HU(H,-),Vti,-)) = J2j M^vW^^dr,, (3.5) 
i=i 
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with the weight 

Wi&r,) = 1 + kiMM + ^{Q(UKO + Q(V)(0). (3.6) 

Here n\ and K2 are two constants to be specified later; Q([/)(£) and Q(V)(£) are the Glimm's total 
wave interaction potentials in U and V at "time" £, respectively; A (£,77) is t ne total strength of 
waves in U and F which approach the i-wave tfe(£,r/) at "time" £, defined by 

A(e,r?) = H i (C,r ? )+C i (e,r ? ) (3.7) 

and 



* = E + E 



a 



^a€j'(!7)Uj'(y) a6j'(!7)Uj'(y)y 
■q a <ri,i<k a <3 r) a >r),l<k a <i 



Ci 



(Y. ae j(U), Va < V ,k a =i + Y. ae j(v), Va > v ,k a =i)\ a \ lf %(??) < 0, 

(Y,a£j(V),r, a < V ,k a =i + T,a£j(U), Va > V ,k a =i)\ a \ lf %(??) > 0, 



where J{U) and J{V) are the sets of fronts at "time" £ in J7 and V, respectively, rj a is the position 
of the front a, and /c a is the characteristic family associated to the front a. In the following, we 
sometimes drop the dependence of £ and write simply qi(r)), Wj (77), Ai(r)), etc. when no confusion 
arises. 

Since, for any £/(£, •)> ^(£> ') ^ BV with TV(C7 (•)) +TV(V(-)) sufficiently small (depending only 
on m, K2, and cf), we have 

3 00 

Co 1 \M£, •) - v& -)|| L i (M+) <E / 1* ^ c ol|i/& •) - -)IUi(R+). 

i=1 

l<Wi(7/)<2, * = 1,2,3, 
for some constant Co depending essentially only on the system and U+. Therefore, for any £ > 0, 

\mz, ■) - v& -)\\ L1{R+) < m v) < •) - v& -)iili(r + ), (3.8) 

where Ci is independent of £, £/", and V. 



3.2. Decreasing of the Lyapunov functional. From now on, we examine how the Lyapunov 
functional ^(£7 5l (£, •), V S2 (£, •)) evolves in the flow direction £ > 0, where U Sl and V^ 2 are approxi- 
mate solutions obtained by the front tracking method from the initial data Uq and Vo, respectively, 
and we also set S = max(<5i,<52) to control the errors. For simplicity, in the following, we also write 
U Sl and V^ 2 as U and V, respectively. 

Denote \i{rj) the speed of the i-wave qi{rj) along the Hugoniot curve in the phase space (it is 
not the characteristic speed except for the characteristic discontinuity). At £ > which is not an 
interaction "time" of the waves either in U or V, for J = J{U) U JiV) (it is a finite set), it holds 
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that 

3 3 

aej i=l i=l 
3 

= E E ( k07a)| Wifaa) (>7a - HVa)) ~ M^)| Wi(j£) (r]a ~ 
a&J i=l 
3 

i=i 

where f] a is the speed of the Hugoniot wave a £ J , b = + stands for the points near the boundary 
r] = 0, and rjb is the slope of the boundary, which is actually zero. For the second equality, we 
used the fact that, when Uq — Vq G if a' is the uppermost front in J, then Xi{rf^,) = 

for i = 1,2,3. Also, suppose that the front /?' € J is the closest to r/ = so that A,(6) = Aj(r/r,). 
Define 

£a,i = k + 1 ^+ (A+ - - | gr I (K ~ Va) , (3-9) 
E b , i = \q i (b)\W i (b)(r ]b + Xi(b)), (3.10) 

where gf = qi (rj±), Wf = W^t), and Xf = A;(?7±). Then 

3 3 

^ # (tf(0, ^(0) = E E E *t + E e h- ( 3 - u ) 

oej i=i j=i 
Our main goal is to establish the following bounds: 

3 

y~]i^6,i < (near the boundary), (3-12) 

i=l 
3 

^^-E'a.i < 0(1)5 \a\ (when a is a weak wave in J). (3.13) 

i=i 

If these are established, from (|3,12[) - (|3.13|) . we conclude 

±$(U(0,V(0)<O{l)5. (3.14) 

If the constant K2 in the Lyapunov functional is chosen large enough, by the Glimm interaction 
estimates, all the weight functions Wi(rj) decrease at each "time" where two fronts of U or two 
fronts of V interact. By the self-similarity of the Riemann solutions, $ decreases at this "time". 
Integrating (I3.14p over the interval [0,£], we obtain 

$ (U(0, V(0) < * (U(0), V(0)) + 0(m (3-15) 

as desired. Actually, (|3.13p has been proved in [3] or [12] via a case by case analysis. 

We also need to consider what happens to the Lyapunov functional $ if a weak wave reflected 
off the boundary. Even in this case, the weights Wi{rj) can be made to decrease across "time" 
£ = r when a 1-wave interacts with the boundary. Using the modified Glimm interaction potential 
Q and Lemma E2] above, it holds that, across "time" £ = r when a weak 3-wave reflected off the 
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characteristic boundary {rj = 0}, both Q and Q decrease, if k + is chosen sufficiently large (see 
Section 2). One sees that the result holds if K2 3> K\ in (|3.6p . What remains to be proved is (|3.12|) 
near the boundary. 

3.3. Estimates near the boundary. In this section, we focus on the estimate of (|3.12p near the 
boundary. It is different from those for the Cauchy problem. We exploit the exclusive property of 
the boundary condition in (|2.17p : The flows U and V have the same pressure near the boundary. 
Then we construct a piecewise constant weak solution only along the Hugoniot curves determined 
by the Riemann date U(b) and V(b), which are the states of U and V near the boundary (that 
is, the point b is some fixed point (^6,%) near the positive £-axis with r/j, > 0). Let hi(b) be the 
strength of the i-th shock in the Riemann problem determined by U(b) and V(b), and Aj be the 
corresponding wave speed. Then, as the second-family is linearly degenerate, we infer that A2 = 0. 
Recalling that fjb = 0, we conclude 

E K2 = 0. (3.16) 
To estimate E b \ + E{, 3 , we consider the following two cases. 

Case 1: h\{b) = 0. If h\(b) = 0, then there is no first-family wave in the Riemann solution to the 
Riemann problem (U(b), V(b)). Since the second-family waves are the characteristic discontinuities, 
the pressure of the middle-state U m keeps unchanged, i.e. p m = p. As the third-family is genuinely 
nonlinear, there must be a jump of pressure across a wave of the third-family. Thus, in this case, 
there must be no wave of the third- family, that is, h 3 (b) = 0. Hence, q\{b) = = q 3 (b). We 
conclude in this case that 

Eb,i = 0, i = l,3. 

Case 2: hi{b) / 0. For this analyzed in Case 1, one concludes that h^(b) / as 

well. Starting from U(b), go along the 1-Hugoniot curve to reach U\, then possibly along the 2- 
characteristic Hugoniot curve to reach U2, and the 3-Hugoniot curve to reach V(b). Furthermore, 
among the first-family and third-family waves, they must be of distinct type, that is, there can be 
a 1-compressive shock (Lax shock) and a 3-decompressive shock (non Lax shock) or vice versa in 
the Riemann solution to keep the pressure unchanged in the approximate solutions U and V. We 
need to show that 

\h 3 (b)\=0(l)\h 1 (b)\. (3.17) 

If this is true, observing that Ai < and A3 > at U+, so that qi(b) = c^hi{b) with the constants 
c" > to be chosen, we have 

E b>1 + E bi3 = c a 1 \h 1 (b)\W 1 (b)\ 1 (b) + c%\h 3 (b)\W 3 (b)\ 3 (b) 

= (&) I ( - c? PTi (6) I A x (6) I + cf O(l) W 3 (&) I A 3 (b) |) 

< |/ il (6)|(-c?|A 1 (6)|+2c^(l)|A 3 (6)|) 

< 0, 



by using the fact that 



A3 (b) 

Ai(6) 



is bounded (depending only on U+) and then choosing cf quite large. 



This completes the proof of (|3.12p . 
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Then what left is to prove (ETPTj) . We know U x = Si(hi)(U b ) and V b = S 3 (h 3 )(U 2 ). For U close 
to U+, since dSl ^^ | a =0 = r\(U) and especially the third argument of t\ is «i(— Aim) / (it is 
nonzero at the reference state i7+), by the inverse function theorem, we infer that |/ii(6)| ^ \p\ —p\, 
and similarly for the third-family, | /13(b)! ^ \p2 —p\- Since P2 = pi, we conclude that |/ii(6)| ^ 
I /13(b) I . Here, for positive quantities a and b, we use a ^ b to mean that there is a positive constant 
Ci depending only on U+ so that C^ 1 b < a < C±b. 

4. Existence of a Semigroup 

Using the existence and uniqueness results established in the earlier sections, we can now establish 
the existence of the semigroup S of solutions generated by the wave-front tracking algorithm. 

Proposition 4.1. Assume that TV (U(-)) is sufficiently small. Then, for 5 > 0, the map (£/(•),£) 1— > 
U s (£,-) := 5|(I7(-)) produced by the wave-front tracking algorithm is a uniformly Lipschitz contin- 
uous semigroup such that 

a) s s u = u, 44^=4+^ 

Proof. Property (i) follows immediately because S s is produced by the wave-front tracking method. 
We now prove property (ii). 

Let II s and V s be the front tracking 5-approximate solutions of (|2.17p with the initial data 
U = Uq(-) and V = Vq(-), respectively. By (|3.8|) and (|3.15|) . we obtain that, for any £ > 0, 



u s (o-v 5 (0 < c$(u 6 (0,v 6 (0) 



T 5 (S),V S ( 



< C$(C/ a (0),F a (0)) + CiC>(l)<5£ 

< C^WV -V\\ L1 + +C l O{l)5i. 



This establishes the Lipschitz continuity of the ^-semigroup with respect to the initial data and 
time. □ 

Definition 4.1. Given 8q > 0, define a domain @ as the closure of the set consisting of the functions 
U : R+ i-> M 3 such that U belongs to L X (R+; M 3 ) by modulo a constant and BV(U-U+)(R + ) < e . 

The semigroup S generated by the wave-front tracking method is given by the subsequent theo- 
rem. 

Theorem 4.1. Let TV(U(-)) be sufficiently small. Then the sequence S s generated from the front 
tracking algorithm is a Cauchy sequence in the L^-norm, and the sequence S$(U) converges to a 
unique limit Sc(U) as 5 — > 0. The map S : [0,oo) x @ 1— >■ Q) is a uniformly continuous semigroup. 
In particular, there exists a constant L such that, for all £1,^2 > and U, V £ 3>, 

(i) Semigroup Property: SqU = U, = S^ 1+ ^ 2 U; 

(ii) Lipschitz Continuity: \\S^U — S^V || L i( R +) < L \\U — ^|| i i( R +)/ 

(iii) Each trajectory £ 1— > S^U yields a weak solution to the initial-boundary value problem (|2.17j) ; 

(iv) Consistency with Riemann Solver: For any piecewise constant initial data U € 
there exists a small 9 > such that, for all £ G [0,0], the trajectory [/(£, •) = S^U agrees 
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with the solution of (|2.17p obtained by piecing together the standard entropy solutions of 
the Riemann problems. 

Theorem 14.11 implies that a sequence of ^-approximate front tracking solutions to (|2.17p converges 
to a unique entropy solution whose value is in @, and this unique limit is L 1 -stable. More precisely, 
we have the following immediate consequence from Theorem 14,11 

Corollary 4.1. Let TV(f7(-)) be sufficiently small. Then the entropy solution to the initial- 
boundary value problem ()2. 17j) produced by the wave-front tracking algorithm is L 1 -stable and 
unique. 

The proof of Theorem 14.11 relies on the subsequent crucial estimate (also used in [3]) about the 
approximation of Lipschitz flows: 

Suppose that S : [0, oo) x <2l — > Q> is a global semigroup with a Lipschitz constant L. Let T > 0, 
Z£f, and Z : [0, T] t— )■ Q be a continuous mapping taking values in piecewise constant functions, 
with jumps along finitely many polygonal lines in the (£, ??)~plane. Then 

\\Z(T)-S T Z\\ Tl <L.{\\Z(0)-Z\\ rl+ F lim 11^ + M) ~ ^(011^ 1 (41) 

Using the essential estimates shown in the earlier sections, Theorem 14, II is shown along a similar 
line of arguments to the one followed by Bressan-Colombo in [3] . The only significant difference here 
between the wave- front tracking by Bressan (cf. [2]) and the version of the method (cf. Holden and 
Risebro |12| ) we employ is on how to control the number of fronts from growing to infinity within 
finite time. In [3], the accurate Riemann solver (ARS) and simplified Riemann solver (SRS) are 
used to construct the approximate solutions. With (ARS), waves of each family can be possibly 
introduced in the Riemann solution where every rarefaction wave is divided into equal parts to 
obtain a rarefaction fan of wave-fronts, while with (SRS), all new waves are lumped together as a 
single non-physical front, traveling faster than all wave speeds. Furthermore, in [3], the simplified 
Riemann solver, rather than the accurate Riemann solver, is employed when the interaction term is 
less than a cut-off function in the order of y/5. In the front tracking method we use, the approximate 
Riemann solution is changed by removing weak fronts. The point of this procedure is that the fronts 
of high generation are quite weak. Particularly, in our arguments, when the interaction term is 
less than 5, this corresponds to new and removed fronts with generation higher than N (here iV 
is a suitably large positive integer computed using the initial approximation parameter 5, see [5] 
and [12J for further details). We also refer to Chen-Li [6] for a related analysis with the full Euler 
system where (SRS) is used in place of (ARS) when the interaction term is less than 5. 

Thus, using an argument similar to [3!|, one concludes that S^"Z n is a Cauchy sequence, converg- 
ing in the L 1 -sense as long as the error from removing weak fronts tends zero as the initial error 
parameter S — > (this has been proved in [5]). Consequently, the map S : [0, oo) X \-± & defined 
as the limit of the approximate solutions produced by the front tracking algorithm is well-defined. 

Now, the proofs of statements (i) to (iv) in Theorem 14. II are as follows. Statements (i), (ii), and 
(iv) are immediate since S is the limit of front tracking approximate solutions S s . It is similar 
to prove (hi) as [3], but, as discussed above, the wave- front tracking method we employ here is 
slightly different. Lastly, we observe that the entropy solution fulfills the boundary condition of the 
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initial-boundary value problem (|2.17p due to the construction of our approximate solutions. This 
completes the proof of Theorem 14.11 

5. Uniqueness in the class of viscosity solutions 

In this section, we analyze an arbitrary Lipschitz semigroup defined on the domain of BV 
functions. In particular, we prove that the semigroup S is uniquely determined when the local flow 
is assigned in connection with the initial data, a piecewise constant function. Based on the results 
in Section 4, we first show that the semigroup S generated by the wave-front tracking method is the 
canonical trajectory of the standard Riemann semigroup (SRS). Then the uniqueness of entropy 
solutions is shown to extend to a broader class, namely the class of viscosity solutions as introduced 
by Bressan in [I]. The essential part here is to show that, in the viscosity class, the entropy solution 
matches the semigroup trajectory generated by the front tracking method. 

Definition 5.1. Problem (|2,17p is said to admit a standard Riemann semigroup if, for some small 
£o > 0, there is a continuous map TZ : [0, oo) x *2l \- > £t and a constant L such that for every 
U,V G *2i and 6, 6 > we have: 

(i) TZ U = U, n^TZ i2 U = K Cl+ s 2 U; 

(ii) ||7^I7-^F|| L1 < L||I7-F|| L1 ; 

(iii) If U £ S 1 is piecewise constant, then, for all £ £ [0, £o] sufficiently small, the trajectory 
U(i, •) = S^U coincides with the solution of (|2.17p obtained by patching together the 
standard entropy solutions of the Riemann problems produced by the discontinuities of U. 

Theorem 5.1. Suppose that the initial-boundary value problem (|2.17p admits a standard Riemann 
semigroup TZ : [0, oo) x Si \— > Si . Consider S the semigroup generated by the front tracking algorithm, 
that is, S(:(U) = lim 5 ^o 5|(tT) with U £ £>. Then, for alii > and U G 9, K(U = S(U. 



Theorem 15.11 is proved using similar arguments as in pQ based on the fundamental estimate (|4.ip 
and the essential feature of the local flow in the ^-direction, that is, the wave-front tracking method 
and the standard Riemann semigroup both have the structure of the Riemann solutions. 

In what follows, we discuss some necessary and sufficient conditions for a function £ i— > U(£) 6 
*3l to coincide with a semigroup trajectory. Following Bressan [I], there are two types of local 
approximate parametrices for system (|2.7p . 

The first approximate parametrice comes from the self-similar solution of the Riemann problem. 
To that end, consider a function U : [0, oo) x M + i-> M 3 and a fixed point (r, Q in the domain of U . 
Suppose that U(t, •) G Q). The boundedness of the total variation implies that the following limits 
exist: 

U~ = lim U(r,rj), U + = lim C/(r,r/). 

Let i? = rj) be the solution of the Riemann problem with piecewise constant data U~ and U + 
and, for £ > r, define the following function: 




C) if |r/-C| < A(e-r), 
if |77-C| > A(e-r). 
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Here A is an upper bound for all wave speeds, 

sup|A fc (f7)| < A, k = 1,2, 3. (5.1) 

U 

The other required parametrice is obtained from the corresponding quasilinear hyperbolic system 

A(U)d ( U + B(U)d v U = 

by "freezing" the coefficients of the matrices A{U) and B{U) in a neighbourhood of the state 
U(r,C). 

For £ > r, define Hnj T Q as the solution of the linear Cauchy problem with constant coefficients: 

Ad^M + Bd n M = 0, M(t, rj) = U(t, rj) . (5.2) 

where A := A(U{t, ()) and B := B{U{t,C)) the matrices evaluated at the fixed state U{t,Q). 
Note that the functions and H depend on the values U (r, () and U (r, C±)- I* 1 the remaining, 
we introduce the class of viscosity solutions and discuss how these solutions share the same local 
characterization as H# and . 

Definition 5.2. Assume that U : [0, T] i— > 3) is a continuous map with respect to the Z^-norm. 
Then the function U is called a viscosity solution of problem (|2.17p if there are constants C and A 
providing bound (|5.ip such that, with f3 and 5 small enough, we have, for every (r, () G [0, T) x M + , 

C+/3-5A 

\ J \U(T + 5,r ] )-H*^ ) (C,v)\^ < CTV{C/(r):(C-/3,C)U(C,C + /3)}, 

C-/3+5A 

J / ^(r + ^r?)-^^)^^)!^ < C(TV{[/(r):(C-/3,C + /3)}) 2 - 

C-I3+5X 

Theorem 5.2. Suppose that problem (|2.17p admits a standard Riemann semigroup 1Z. Then a 
continuous mapping U : [0, T] i— )■ St is said to be a viscosity solution of (|2.17p if and only if 

U(£,-) = Tl^J at every £ G [0,T]. (5.3) 

Using Theorem 5.2, one concludes that, for problem ([2.17p . the entropy solution is unique in the 
viscosity class of solutions, which coincides with the semigroup trajectory S^U generated by the 
front tracking method. More precisely, we have 

Corollary 5.1. A continuous mapping U : [0, T] i— > *2J is a viscosity solution if and only if 

U(£,-)=S ( :U for any £ G [0,T]. (5.4) 

With the estimates in Sections 2-3, the proof here follows along a similar line of reasoning to the 
one presented in PQ. The only difference is a straight boundary and physical domain restricted to 
the positive 77-axis; nonetheless, we can still proceed with the proof provided that the convergence 
of the wave-front tracking method is obtained which has been outlined in Section 2 and carried out 
in Eulerian coordinates in [5]. 
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Remark 5.1. Note that, in simpler cases such as the potential flow, isentropic or isothermal Euler 
flow, as far as the L 1 -stability problem is concerned, we achieve the same results as the full Euler 
equations. 

6. Uniqueness of Solutions to the Free Boundary Problem in Eulerian 

Coordinates 

In this section we apply the following Wagner's Theorem |X5|, Theorem 2] to show the uniqueness 
of entropy solution to the free boundary problem (|1.7|) . 

Theorem 6.1 (Wagner [IS]). Consider 

d t U + d x F(U) = 0, (6.1) 

where (t,x) £ M 2 , U(t,x) = («!,-•• ,u n ) T , and F(U) = (fi{U),--- , f n (U)) T 6 R n is a smooth 
vector-function. For any bounded measurable solutions of (|6.ip . with u±(t,x) > 0, let (t,y) satisfy 

%± = u 1 (t,x), % = -h{U{t,x)) 
ox ot 

in the sense of distributions. ThenT : (t, x) \— > (t,y(t,x)) is a Lip schitz- continuous transformation, 
which induces a one-to-one correspondence between L°° weak solutions of (16. ip on R + xl satisfying 
< e < u\(t, x) < M < oo for some e and M , and L°° weak solutions of 

Ul Ul 

dt{—(u 2 ,--- ,u n ) T )+d y {(f 2 (U),--- J n (U)) T -— fx(U)(u 2 ,--- ,u n ) T ) =0, (6.2) 

Ul Ul 

on M + x R satisfying e < u\(t,x) < M. 

If rj(U) is any convex extension of (|6.ip . i.e., there is a flux q(U) such that Vr/VF = X7q, so 
that dtr] + d x q = for classical solutions, then any solution of (16.11) . satisfying 

d t r](U) + d x q(U) < (6.3) 

in the sense of distributions in (t,x) corresponds to a solution of (16. 2p . satisfying 

d t fj(V) + d y q(V) < 0, (6.4) 

in the sense of distributions in (t,y), where V = (vi, ■ ■ ■ ,v n ) T , fj(V) = ~p-, and q{V) = q{U) — 
fi(U)fj(V) . Furthermore, rj is convex if and only iffj is convex as a function ofV. Thus, the Lax's 
entropy inequality holds for a solution of (|6.1|) if and only if it holds for the corresponding solution 

of ([52]). 

As indicated by Wagner |15t p. 123], this theorem may also be applied to initial-boundary value 
problems by slight modification. Thus, we can also use it in our case. Suppose that Uq(d) and 
Uq{v) are two initial data functions of problem (|1.7p in Eulerian coordinates, and {U l ,g l ) and 
(U 2 ,g 2 ) are the corresponding weak entropy solutions obtained by the front tracking method. We 
have known that U 1 and U 2 must be bounded. We now show that U 1 = U 2 and g 1 = g 2 , provided 
TJ 1 - IT 2 

Suppose that V 1 and V 2 are the solutions of (|2.17p in Lagrangian coordinates that correspond 
to U 1 and U 2 , respectively. By Theorem 16.11 it suffices to show that V 1 = V 2 . Then, from (|2. 18j) . 
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we may find g 1 = g 2 . Therefore, by the uniqueness results of problem (|2.17p (Corollary 14. ip . it 
suffices to prove the following lemma. 

Lemma 6.1. IfU&(y) = U${y), then V l {Q,r]) = V 2 (0,7]). 

Proof. For j = 1,2, let T 3 : (x,y) t— >■ (£,77) be the Lagrangian transform (|2.4p associated with the 
solution U 3 . They are Lipschitz-continuous and one-to-one. Furthermore, we know U J (x,-) is of 
bounded variation and lim x -+ Xo CP (a;, •) — U 3 (xq, -)|| ri = 0. Thus, we can not only solve r] from 
(|2.2p . but also those equations make sense on each line x = constant. 
We first show that, if U$(y) = U%(y), then T 1 ]^ = T 2 \ x=0 . 

Suppose (£j,T)j) = T J (x,y). Recall by definition that f = x. Then, at x = 0, we have both 
£i = £2 = 0. Note that rjj satisfies 7^(0,0) = and solves 



x =o = (pv,y(0,y). 



dy 

Since (pu) 1 (Q, y) = (pu) 2 (0,y) from U^y) = U^y), we conclude that 771(0,2/) = 772(0,7/) as well. 
This shows that T 1 (0,y) = T 2 (0,y), which implies that (T 1 ) _1 (0,?7) = (T 2 ) _1 (0, 77) and lies on 
{x = 0}. Since 

V j (t;,v) = Ui((Tir 1 (t;,r } )), 

we conclude 



V\0, V ) = UW 1 )- 1 ^)) = UHiT 2 )- 1 ^^)) = U 2 ((T 2 )-\0, V )) = F 2 (0,77). 



Finally, summarizing all the analysis above, we state the main theorem of this paper. 



□ 



Theorem 6.2. For the initial data Uq close to the reference state V '+ in the sense that ||[/q — C/+|| BV 
is sufficiently small, there is one and only one weak entropy solution (U,g) to problem (|1 .7j) con- 
structed by the front tracking method. Furthermore, reformulating this problem in Lagrangian 
coordinates as problem (|2,17p . then the L 1 -stability holds in the sense that 

\\v\0 - ^ 2 (£)|| L1{R+) < c \\v\o) - v 2 (o)\\ L1(R+) 

for any two solutions V 1 and V 2 of (j2. 17j) . The solution is also unique in the class of viscosity 
solutions in Lagrangian coordinates. 
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